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Summary. In the framework of ion-trap quantum computing, we de- 
velop a characterization of experimentally realistic imperfections which 
may affect the Cirac-Zoller implementation of the CNOT gate. 
The CNOT operation is performed by applying a protocol of five laser 
pulses of appropriate frequency and polarization. The laser-pulse proto- 
col exploits auxiliary levels, and its imperfect implementation leads to 
unitary as well as non-unitary errors affecting the CNOT operation. 
We provide a characterization of such imperfections, which are physi- 
cally realistic and have never been considered before to the best of our 
knowledge. Our characterization shows that imperfect laser pulses un- 
avoidably cause a leak of information from the states which alone should 
be transformed by the ideal gate, into the ancillary states exploited by 
the experimental implementation. 

Key words: Ion-trap quantum computing, Cirac-Zoller CNOT gate, decoerence, 
error models. 



1 Decoherence and Quantum Computing by Ion Traps 

Quantum computation and quantum communication require special physical 
environments: Decoherence, noise and experimental imperfections threaten the 
correctness of quantum computations, as well as the intended functioning of 
quantum communication protocols. Several experimental implementations have 
been proposed and explored in order to build systems capable of having limited 
sensitivity to unwanted perturbations, while allowing both desired interactions 
among internal components and access from external systems or users. 

Ion-trap quantum computation is one of the main and rapidly evolving ex- 
perimental possibilities. Several new ideas and experimental techniques [1] [2] [3] 
suggest that ion traps offer a promising architecture for quantum information 
processing. In this framework, we develop a characterization of experimentally 
realistic imperfections which may affect the Cirac-Zoller implementation of the 
CNOT gate [3]. The CNOT operation is performed by applying a protocol of five 
laser pulses of appropriate frequency and polarization. The laser-pulse proto- 
col exploits auxiliary levels and this could origin unitary as well as non-unitary 
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errors affecting the CNOT operation. We provide a characterization of such im- 
perfections, which are physically realistic and have never been considered before 
to the best of our knowledge. 

The paper is organized as follows. In Section[2J we briefly review the sequence 
of five laser pulses which implement the CNOT gate by ion-trap techniques, as 
proposed by Cirac and Zoller. In Section [3l we explore the physically realistic 
perturbations which may affect the laser-pulse protocol for the Cirac-Zoller im- 
plementation of the CMINUS gate and, consequently, of the CNOT gate. First, in 
Section 13.11 we describe by diagrams of quantum states the twelve-level Hilbert 
space in which we reproduce both the ideal and imperfect actions of the CNOT 
gate. Subsequently, in Section [3~T2l we model imperfections in all the three laser 
pulses which implement the CMINUS gate, by applying perturbations on the pa- 
rameters determining the physical characteristics of the lasers. Then, in Section 
I3.3l we formulate our error model by means of the well-known density-matrix for- 
malism and Kraus representation, showing how leakage errors are unavoidable 
in realistic conditions. Finally, in Section [4] we express our conclusive remarks. 

Throughout this paper, we assume the reader is familiar with the basic ion- 
trap techniques. Useful and detailed physical and computational descriptions of 
ion-trap quantum computing can be found, for instance, in [5] (pages 544-561) 
or 0. 



2 The Cirac-Zoller CNOT Gate 

Ion-trap quantum gates can be obtained by applying to ions implementing qubits 
the appropriate combination of laser pulses, tuned to the appropriate duration 
and phase. In controlled two-qubit operations, the motion state of the string of 
ions (the phonon qubit) is exploited as a "bus" to transfer quantum information 
between two qubits implemented by ions. 

Following the proposal of Cirac and Zoller [4] , we briefly recall how to obtain 
a CNOT quantum gate (CNOTcz) acting on a set of N ions, with the ion / as the 
control qubit and the ion m as the target qubit. 

A general state in the computation is described as \ii, i m ; i n ), where the 
last position is always filled by the state of the phonon qubit; for the sake of 
simplicity, we omit all the qubits of the register which are not affected by the 
protocol. Qubits can be in the ground or excited levels \g) and |e); at some point 
of the protocol, an auxiliary level |o) is also necessary. The initial state for the 
computation is \ii,i m ;0). 

First, a red de-tuned laser acts on the ion I in order to map the quantum 
information of the control ion onto the vibrational mode. This impulse changes 
the states of the two-qubit computational basis as follows: 

\9h9m;0) -> \gi,g m ;0) 

|# ; ,e m ;0) -> \gi,e m ;0) ,^ 
\ei,g m ;0) -> i\gi,g m \ 1) 
|e;,e m ;0) -> i\gi,e m ; 1). 
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Then, a red de-tuned laser is applied to the ion m. The impulse is now expressed 
in the basis {\g m ]n — 1), \a m ;n — 0)}, thus exploiting the auxiliary level \a m ), 
and it acts as follows: 



' \gi,g m ;0) 

i\ei,g m ; 1) 
, i\gi, e m ; l) 



\gi,g m ;0) 
\gi,e m ;0) 

-i\ei,g m ; 1) 
i\9h e m] l)- 



(2) 



Finally, a red de-tuned laser acts once again on the ion I, in order to map the 
quantum information of vibrational mode back onto the control ion: 



\gi,g m ;ty 

\gi,e m ;0) 
-i\gi,g m ; l) 
i\gu e m; i) 



\gi,g m ;0) 

\gi,e m ;0) 
\ei,g m ; 0) 

-|e;,e m ;0). 



(3) 



The global effect of the three laser pulses is to induce a controlled phase- 
shift gate with a phase-shift of an angle S = it, that is, a CMINUS gate. From 
the CMINUS gate, the CNOT gate can be obtained by applying a single-qubit 
Hadamard gate before and after the three laser pulses. 

In conclusion, the CNOTcz operation between the target qubit m and the 
control qubit I can be performed by applying a protocol of five laser pulses of 
appropriate frequency and polarization. 



3 Imperfections in the Cirac-Zoller cnot Gate 

In the ion-trap implementation of the CNOT gate proposed by Cirac and Zoller, 
a key-role is played by the laser protocol implementing the CMINUS operation, 
which exploits ancillary levels. The role of such additional levels is here explored 
to model unitary as well as non-unitary errors affecting the CMINUS operation 
and, consequently, the whole implementation of the CNOT gate. 

The error model schematically consists in the following steps. First, unitary 
errors are applied to each impulse gate constituting the CMINUS protocol. Then, 
a partial trace operation is performed on the collective vibrational motion, whose 
levels are neglected at the end of both ideal and non-ideal computations. In the 
ideal protocol, a final step of partial tracing on the ancillary levels would lead to 
perfectly reproducing the action of the CMINUS operation. However, in a general 
experimental situation, our error model shows that it is no longer possible to 
trace over the ancillary level without loss of information: The unitary errors 
introduced in the three CMINUS laser pulses unavoidably cause an irreversible 
leak of information in the ancillary states. 

3.1 The Hilbert Space 

In order to describe the most general perturbation of the CMINUS laser-pulse 
protocol for the Cirac-Zoller implementation of the CNOT gate, we rely onto a 
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useful representation and ordering of the twelve-level space which reproduces 
the action of the gate. 

This ordering is illustrated by the diagram of states in Figure [T] In the usual 
way of representing quantum computations by quantum circuits, each horizontal 
line represents a qubit. In diagrams of states [7], we draw instead a horizontal 
line for each state of the computational basis, here adding horizontal lines for the 
additional levels necessary to reproduce the action of the CNOT gate. Similarly 
to quantum circuits, any sequence of operations in diagrams of states must be 
read from left (input) to right (output). 

We first represent the ideal action of the three laser pulses, that is, we repre- 
sent how the information contained in the states of the considered twelve-level 
space is elaborated when the CNOT gate is implemented by a system unaffected 
by errors or noise. We put the control qubit m, which has two possible levels 
and 1, in the most significant position; we put the target qubit I, which has 
two possible standard levels and 1 plus an additional ancillary level 2, in the 
middle; finally, we put the phonon qubit n, which has two possible levels and 
1, in the least significant position (see upper scheme in Figure [T]). We order the 
resulting twelve levels accordingly (see lower diagram of Figured]). The ordering 
is such that in the top half of the diagram we have the six levels corresponding 
to the non-excited state of the phonon qubit, namely, n = 0; in the bottom half 
of the diagram we have the remaining six levels, corresponding to the excited 
state of the phonon qubit, namely, n — 1. Finally, the two bottom levels of each 
half of the diagram correspond to the ancillary level of the target qubit, namely, 
m = 2. 

Since we are considering initial states for which i) the collective vibrational 
motion is not excited (n = 0) and ii) the auxiliary level stores no information, 
the initial information is stored only in the first four states |000), |010), 1 1 00) 
and 1 110). The adopted ordering allows us to collect at the top of the diagram 
the four levels containing the input information, when the collective vibrational 
motion is not excited and the ancillary states have not yet been exploited. After 
applying in sequence the three operations corresponding to the three laser pulses 
previously described, we show at the rightmost end of the diagram how the global 
transformation affects the twelve levels of the overall system. Again, at the top 
of the diagram we read on the first four states how the computation reproduces 
the action of an ideal CMINUS gate. 



3.2 Imperfect Impulse Gates 

Referring to the Hilbert space previously described, we now model imperfections 
in the three impulse gates involved in the implementation of the CNOTcz gate. 
The action of a general unitary matrix acting on a two-qubit system, 

jj _ cos I (cos -0 + i sin ip) i sin | (cos + i sin (j>) 
i sin | (cos (j> — i sin 4>) cos | (cos ip — i sin i[>) 

is determined by the three parameters 9, ip, <fi. When considering unitary trans- 
formations induced by laser pulses, these parameters are associated with precise 
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|n) - phonon qubit (2 levels) 

|m) - target qubit (3 levels) 

|1) - control qubit (2 levels) 




Fig. 1. Diagram of states [7] representing a noisy ion-trap implementation of the 
CMINUS gate, according to the Cirac-Zoller protocol. The control qubit m in set the 
most significant position, with possible levels and 1; the target qubit I in set the 
middle, with possible levels 0, 1 and 2 (ancillary); the phonon qubit n in set the least 
significant position, with possible levels and f . The ordering of the twelve resulting 
levels is such that the top six levels of the diagram correspond to the non-excited state 
of the phonon qubit (n = 0), while the bottom six levels of the diagram correspond 
to the excited state of the phonon qubit (n = 1). In each half of the diagram, the 
two bottom levels correspond to the ancillary state of the target qubit (m = 2). The 
sequence of operations must be read from left (input) to right (output); information 
flows on the thick lines, while thinner lines correspond to absence of information. Note 
that the action of the CMINUS gate can here be visualized by simply following the thick 
lines corresponding to the first four levels, which alone store all the input information 
at the start of the gate operation. 



physical features of the impulses: The parameter 6 — fit denotes the impulse 
area, where fl is the intensity and t is the duration of the impulse; the parameter 
ip is related to the laser de-tuning; the parameter (f> denotes the laser phase. By 
appropriately varying these parameters, we obtain the unitary matrices which 
correspond to the ideal impulses implementing the CMINUS gate. 

The first and the third impulse matrices are obtained by imposing: 



{6 = 7T,V = 7T, = 0}, C/ (1 < 3) = 



i 

1 



(5) 
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while the second impulse matrix is obtained by imposing: 



{6» = 7r,V' = 7r,0 = O}, U {2) = 



-1 
-1 



(6) 



We apply a perturbation on all the three parameters of all the three impulses, 
which we will denote from now on with index r. 



{6' + A8\ ift + Zty/, ft + Aft} 
Thus, each noisy impulse is now described by: 

[/« = 



i = 1,2,3. 



u \\ u 12 
u 2l u 22 



(7) 



(8) 



where, for i = 1, 2, 3, we have: 
L + A0 L 



It, i = cos ■ 



(cos(t/>' + Aip L ) + ism(ip L + Aip 1 )), 



u\ 2 = i sin (cos(0 t + Aft) + isin(ft + Aft)), 



u 2 i = i sin ■ 



2 

n L + A9 l 
~2 
A0 l 



(cos{ft + Aft) - i sm(ft + Aft)), 



U 2 2 = COS ■ 



(cos^ 1 + Aip L ) -ism(ip L + Aip L )). 



(9) 



All the three noisy impulse matrices are subsequently embedded in the 
twelve-level space of the Cirac-Zoller implementation, obtaining respectively the 
matrices F«, and The entries {u\^} are given by the parameters 

defined by equation ([5]); since the matrices are very sparse, we omit all their null 
entries to obtain a clearer visualization of each matrix structure: 



1 . 
1 



1 



'ii 



'12 



'■21 



'21 



'22 



'22 



VP) 



'12 



1/(3) = 
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1 



"11 



'21 



'21 



t»12 



'22 



'22 



1 



(10) 



Finally, we appropriately embed the Hadamard gates in the twelve-level space 
of the Cirac-Zoller implementation (still denoting by H the resulting matrices 
of dimension 12 x 12, for the sake of simplicity). 

By multiplying the five matrices in the appropriate order, we obtain the 
overall matrix of dimension 12 x 12 which describes the action of the CNOTcz 
gate: 

CNOTcz = H V {3) V {2) V (1) H. (11) 



3.3 Density Matrix Evolution and Kraus Operators 

We now calculate the density matrix transformation corresponding to the pre- 
viously described model for imperfections in the impulse gates originating the 
CNOTcz gate. 

The initial information is stored only in the first four states 1 000) , |010), 
1 100) and 1 110) : This is equivalent to say that all initial information is stored in 
a density matrix of dimension 4x4, which we denote by p l ^ x ^. Thus, the overall 
initial density matrix PT2X12 nas non-zero entries only in the 4x4 upper-and- 
lcftmost positions: 



P 6x6 — 



P4x4 u 4x2 
2x4 2X 2 



P 12x12 



6x6 6x6 
06x6 0g x6 



(12) 



The evolution of the overall initial density matrix caused by the laser-pulse 
protocol is given by: 



P 12X12 = CNOTcz P T 2x 12 CNOT cz . 



(13) 



In order to determine the Kraus operators describing the evolution of the reduced 
density matrix, the unitary matrix CNOTcz can be decomposed into sub-matrices 
A, B, C and D, each one of dimension 6x6. Thus, equation (fl~3| can be written 
as follows: 



fin 
P 12x12 



A B 
C D 



P 6x6 
6x6 



6x6 
6x6 



At ct 



(14) 



Since at the end of the process implementing the CNOTcz gate we neglect 
the levels corresponding to the collective vibrational motion, we now perform a 
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partial trace operation on the overall final density matrix p 12x12: m respect to 
the phonon qubit n: 



P 6x6 — Tr n = LSB { p 12x12 } — ^ p'gxg A* + Cp'g x6 C^. 



(15) 



We further decompose the sub-matrices A and C, both having dimension 6x6, 
into sub-matrices {Ai,Ci}, for i = 1 . . .4, each one having the appropriate di- 
mension as expressed in the following equation: 



.fin 
P 6x6 



Al 4x4 A2 4x2 
^43 2x4 A4 2x2 



P4x4 4x2 
2x4 2x2 



41 T 4r 

4x4 ^° 4x2 



A2 



2x4 



A4 



2x2 



+ 



+ 



CI 4x4 C2 4x2 
C3 2x4 C4 2x2 



P4x4 4x2 
2x4 2x2 



oi 4x4 °°4x2 
2x4 C4 2X 2 



(16) 



Consequently, the final reduced density matrix can be expressed as: 



fin 
P 6x6 



4x4 

p2 

4x2 

p3 

2x4 2x2 



(17) 



where: 



Pi 4x4 — ^1 4x4 ^'4x4 ^1 4x4 + CI 4x4 ^'4x4 4x 4, 



M 2x2 — A3 2x4 ^4x4 ^3 4x2 + C3 2x4 P4x4 C3 4x2 



(18) 
(19) 



Recalling expression (|11[) . we finally explicit the Kraus operators j41, ^43, CI 
and C3: 



Al 



4x4 



10 
10 

00o|a^ 
a\ a\ 



where: 



.43 



2x4 








(20) 



(2«i 1 «? 1 +t4i(l + tiii)«i 2 ), a\ 



4i (- 1 



3 

12- 



(21) 



with parameters {mJ^}, for t = 1,2, 3, defined by equation ([5]), and: 



CI 



1x4 



00c} c\ 
QQc\ -c\ 





C3 



2x4 



c 3 c 3 




(22) 



where: 



c i = -7= ( u ii u ; 



U 21 M ll U 22)i 



V2 

with parameters {wj ■}, for 



2 — yTT ( u ll u 21 + u 21 u 22)7 



1 



Uo i tin 



V2^ 21 21 ' 
1, 2, 3, once again defined by equation ©. 



(23) 
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In conclusion, when the three laser pulses are ideally implemented without 
perturbation, the overall transformation coincides with the ideal CNOT gate, 
that is, the upper-and-leftmost final density sub-matrix p ma i n = /0I4X4 is the 
result of a perfect CNOT gate application to the density matrix of the four initial 
states. On the other hand, when we deal with inaccurate laser pulses, the final 
information is no longer stored only in the evolved four initial states. Errors 
unavoidably cause a spreading of information, which partially remains stored in 
the four evolved initial states, and partially leaks in the ancillary levels involved 
in the action of the second laser pulse. 

More precisely, the final density matrix previously calculated shows how in- 
formation is spread by the sequence of imperfect operations. The density sub- 
matrix Pmain still holds the information corresponding to the imperfect action of 
the gate on the desired control and target qubits, while the lower-and-rightmost 
final density sub-matrix p anc = p4 ^"2 holds the information leaked in the aux- 
iliary levels and thus lost from the control and target qubits. 

The traces of sub-matrices /O ma in and p anc give the probabilities p ma m and 
Pane that information spreads and populates the main and ancillary levels, re- 
spectively. Obviously, p ma i n +p an c = 1, as the overall density matrix has unit 
trace. 

4 Conclusive Remarks 

We have developed a characterization of experimentally realistic imperfections 
which may affect the Cirac-Zoller implementation of the CNOT gate. 

In the framework of ion-trap quantum computing, the CNOT operation can 
be performed by applying a protocol of five laser pulses of appropriate frequency 
and polarization. The laser-pulse protocol exploits auxiliary levels, and this may 
results in unitary as well as non-unitary errors affecting the CNOT operation. 
We have provided a characterization of such imperfections, which are physically 
realistic and have never been considered before to the best of our knowledge. 
This characterization shows that it is no longer possible to disregard the ancillary 
levels without loss of information: Unitary errors introduced in the laser pulses 
unavoidably cause a leak of information in the ancillary states, which can no 
longer be completely retrieved at the end of the computation. 

Modeling experimentally realistic perturbations affecting two-qubit universal 
operations allow us to explore the impact of two-qubit imperfections on quan- 
tum computation and communication protocols. Our previous characterization 
of general single-qubit errors and its application to a quantum privacy amplifi- 
cation protocol based on entanglement purification [8] proved that different im- 
perfections can affect a quantum protocol very differently. In a similar way, the 
modeled two-qubit imperfections may have unexpected effects on the quantum 
computation and communication protocols simulated in their presence. Exten- 
sive simulations will highlight which error parameters are most dangerous and 
which kinds of perturbations should be expected, especially regarding errors on 
the ancillary levels. 
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Finally, by addressing general two-qubit errors affecting universal two-qubit 
operations, we also aim to extend our previously developed general single-qubit 
error model. Since multi-qubit unitary operations can always be computed by 
equivalent quantum circuits composed of single-qubit and two-qubit gates, a 
compact model comprising single-qubit and two-qubit errors would offer a sat- 
isfiable characterization of general imperfections. Unfortunately, describing and 
understanding all parameters which characterize the most general transforma- 
tion of a two-qubit density matrix can not be achieved by straightforward ana- 
lytical study. Thus, exploring imperfect implementations of universal two-qubit 
operations in the main current experimental frameworks can offer a valuable 
help to determine the role and the physical meaning of significant parameters 
describing the most general two-qubit noise transformation. 
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